Spin relaxation in a GaAs quantum dot embedded inside a suspended phonon cavity 
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The phonon-induced spin relaxation in a two-dimensional quantum dot embedded inside a semi- 
conductor slab is investigated theoretically. An enhanced relaxation rate is found due to the phonon 
van Hove singularities. Oppositely, a vanishing deformation potential may also result in a suppres- 
sion of the spin relaxation rate. For larger quantum dots, the interplay between the spin orbit 
interaction and Zeeman levels causes the suppression of the relaxation at several points. Further- 
more, a crossover from confined to bulk-like systems is obtained by varying the width of the slab. 

PACS numbers: 73.21. La, 71.70.Ej, 63.20.Dj, 72.25.Rb 
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Spin properties in nanostructures have become a 
field of intense research ranging from spin field-effect 
transistor^ spin-polarized p-n junctions^ up to quantum 
spin computersi^ The quantum dot (QD) may be a good 
choice for quantum electronics due to its zero dimension- 
ality, quantized energy levels, and long coherence times 
of spin states.^''' For example, the spin of an electron con- 
fined to a QD can form a qubit.'*'^ However, some scat- 
tering processes will cause the change of the spin states. 
One important process is related to the phonon-induced 
spin-flip resulting from the spin-orbit interaction. This 
affects the time of spin purity in the QD. In order to 
keep the information unchanged, a long relaxation time 
is required. 

In general, the spin-orbit (SO) coupling, which is one 
of the main causes of spin relaxation, is a relevant in- 
trinsic interaction in nonmagnetic semiconductors. It is 
known that there are two different types of spin-orbit 
coupling as QDs are fabricated within semiconductors of 
a zincblende structure. The first one is the Dresselhaus 
interaction, which is due to the bulk inversion asymme- 
try of the latticeiS*2iiS The second is the Rashba interac- 
tion caused by the structure inversion asymmetry. The 
spin-orbit couplings mix the spin states with different 
orientations in the Zeeman sublevelsiSiiii^ii^ and there- 
fore make spin relaxation possible in the presence of the 
electron-phonon interaction. 

Relaxation times of electron spins in a QD have been 
measured by electrical pump-probe experimentsji^ The 
triplet-to-singlet transition with emission of phonons was 
found with corresponding spin relaxation times of about 
200 /iS. Recently, the spin relaxation time in a one- 
electron GaAs QD was measured by a similar electrical 
pump-probe technique j"'^^'^^ - As the magnetic field was ap- 
plied parallel to the two-dimensional electron gas, the 
Zeeman splitting of QD was observed in dc transport 
spectroscopy. By monitoring the relaxation of the spin, 
the relaxation time was found to have a lower bound of 
50 /iS at an in-plane field of 7.5 Tiii 

On the theoretical side, spin relaxation between two 
spin-mixed states in semiconductor QDs has been stud- 
ied recently. However, to the best of our knowledge, 
all previous studies of spin relaxation have concentrated 



on QDs embedded in the bulk material,i2i2aiaiiS2iSi24 

whereas studies of spin relaxation induced by confined 
phonons are still lacking. In this work, we therefore con- 
sider a single QD embedded into a free-standing structure 
(semiconductor slab), where the relevant characteristic is 
the two-dimensional phonon wavevector for the acoustic- 
phonon spectrum as shown in Fig. i |^S42SiSL22i2S Since 
the reduced dimension will enhance the deformation po- 
tential, we will mainly focus on the spin relaxation rate 
induced by the deformation potentialiSSiSSiSi In this pa- 
per we describe the model with spin-orbit coupling. En- 
ergy spectra of the QD can be solved by using an exact 
diagonalization method. We then apply the Fermi golden 
rule to calculate spin relaxation rates for typical param- 
eters. We discuss the dependence of the spin relaxation 
rates on the size of the QD, the phonon bath tempera- 
ture, and the width of the slab. 
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FIG. 1: (Color online) (a) Schematic view of single QD em- 
bedded in the semiconductor slab with a width of a. (b) The 
side view shows a QD is located at z = 0. 

We consider an isotropic QD with an in-plane parabolic 
lateral confinement potential. An external magnetic field 
B is applied perpendicularly to the surface of the QD as 
shown in Fig. 1(a). The electronic Hamiltonian of this 
system can be written as 



He — Hq + Hg 



(1) 
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The first term describes tlie electron Hamiltonian with- 
out the spin-orbit coupling, 
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where P = —ihV + (e/c)A is the kinetic momentum 
with vector potential A = (i?/2)(— y, x, 0) confined to 
the 2D plane. Here m* is the effective electron mass, e 
is the electron charge, c is the velocity of light, loq is the 
characteristic confined frequency, g* is the bulk gi-factor, 
/Is is the Bohr magneton, and Uz is a Pauli matrix. 

The Rashba and Dresselhaus interactions 
{Hso = Hr + Hd) are given by 
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'-{-CFxPx + CTyPy). 
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The coupling constants \r and Ad determine the spin- 
orbit strengths, which depend on the band-structure pa- 
rameters of the material. Besides, the Rashba and Dres- 
selhaus terms are also associated to the perpendicular 
confinement field and the confinement width in the z- 
direction, respectively. 

For the electron Hamiltonian Hq, the well-known 
Fock-Darwin states ^n,i,cT can be easily obtained. 
The corresponding electron energy levels are En,i,a = 
Ml{2n+\l\ + 1) + hwBl/2 + o-EsMere n (= 0,1^2...) 
and I (= 0, ±1, ±2...) arc the quantum numbers. The 
renormalized frequency is f2 = \/uJq + w^, with the cy- 
clotron frequency ujb = eB/m* and the characteristic 
confinement frequency uq limited by the effective QD 
lateral length Iq = y/h/m*LOo. Here, Eb = giJ.BB/2 is 
the Zeeman splitting energy, and u = ±1 refers to the 
electron-spin polarization along the z axis. To solve the 
Schrodinger equation with {He = Hq + Hgo), the (spin 
mixing) wave function is expressed in terms of a series of 
eigenf unctions: 'i'l (r, 0) = Cn,l,cr^n,l,cr for each state 
£. After exactly diagonalizing the electron Hamiltonian, 
the corresponding eigenvalues and the coefficient Cn,i,a 
can be obtained numerically. 

Before calculating the spin relaxation rate, the con- 
fined phonon in the free-standing structure must be in- 
troduced here. Following Re. [25], we consider an infinite 
film with width a (Fig. 1). For the effect of the con- 
tact with the semiconductor substrate, we neglect the 
distortion of the acoustic vibrations. Under this con- 
sideration, one can ensure that the in-plane wavelength 
can be shorter than the characteristic in-plane size of the 
solid slab. For simplicity, the elastic properties of the 
slab are isotropic. Small clastic vibrations of a solid slab 
can then be defined by a vector of relative displacement 
u{r,t). Under the isotropic elastic continuum approxi- 
mation, the displacement field u obeys the equation 
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= c2v2u+(cf-c2) V(V-u), 



(5) 



where q and Cf are the velocities of longitudinal and 
transverse bulk acoustic waves. To define a system of 
confined modes, Eq. (5) should be complemented by the 
boundary conditions at the slab surface z = ±a/2. Be- 
cause of the confinement, phonons will be quantized in 
subbands. For each in-plane component q|| of the in- 
plane wave vector there are infinitely many subbands. 
Since two types of velocities of sound exist in the elas- 
tic medium, there are also two transversal wavevectors qi 
and Qt- In the following, we consider the deformation po- 
tential only. This means there are two confined acoustic 
modes: dilatational waves and fiexural waves contribute, 
but shear waves are neglected because of their vanishing 
interaction with the electrons for spin relax;ation. 

For dilatational waves, the parameters qi^n and qt,n can 
be determined from the Rayleigh-Lamb equation 



tan(gt,na/2) _ _^\\%nqt^ 
tan(gi,„a/2) (gji - 
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with the dispersion relation 
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(6) 
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where u>n,qii is the frequency of the dilatational wave in 
mode (n, q||). For the antisymmetric flexual waves, the 
solutions 5; „ and qt,n also can be determined by solving 
the equation 



tan {qi^nd/ 2) MwqLnQt, 



tan (gt,„a/2) {qf - qf j 
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together with the dispersion relation, Eq. (7). 

The electron-phonon interaction through the defor- 
mation is given by H^p = -Badivu, where Ea is the 
deformation-potential coupling constant. The Hamilto- 
nian can be written as 

Hep = ^ MA(q||,n,2;)(a+ -Faq||)exp(iq|| Til), (9) 

q||,n 
X=dJ 

where ry is the coordinate vector in the x-y plane and 
the functions and Mf describe the intensity of the 
electron interactions with the dilatational and fiexural 
waves, and are given by 



Md (gy ,n,z) = Fa,r. 
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X sm(— — )cos(gi,„z) 



(10) 



{q\\,n,z) = Ff^. 
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where A is the area of the slab, p is the mass density, and 
Fd,n {Pf,n) is the the normahzation constants of the rt-th 
eigenmode for the dilatational (flexural) waves. Although 
the fluctuation of the dot (due to strain etc.) may affect 
the spin-orbit and electron-phonon coupling, we, for sim- 
plicity, neglect the effect on the scattering rate in this 
work. 

We calculate the spin relaxation rates between the two 
lowest (spin mixing) states from the Fermi golden rule^ 

r = Y ^ |MAn(/|e^''"-"-"|z)|' 

\=dj 

x{Ng^^+l)6{AE-hiJn^g^^), (12) 

where the energy AE {— Ei — E f) is the energy differ- 
ence between the first excited \i) and ground |/) states. 
A^gii represents the Bose distribution of the phonon at 
temperature T. For the sake of simplicity, we consider 
the QD to be located at z = so that the function Mf 
for flexural waves plays no role. 

400, n 1 



plot the phonon group velocity (dujn,q^^ / dq\\) as a func- 
tion of q|| around the van Hove singularity as shown in 
Fig. 2(b). There are three modes contributing to the 
relaxation rate. In particular, a crossover from positive 
to negative group velocity is observed for one mode. Be- 
cause of the zero phonon group velocity, the rate behaves 
sharply at that magnetic field. However in a real system 
the van Hove singularity would be cut off or broadened 
because of the finite phonon lifetime. Contrary to the en- 
hanced rate, we find a suppression of the spin relaxation 
rate (arrow 2) at small magnetic field (also seen in the 
lower inset). This comes from a vanishing divergence of 
the displacement field u. As can be seen from Eq. (10) in 
detail, the deformation potential disappears at the con- 
dition of — qt (Fig. 2(c)), which causes a zero spin 
relaxation rate. Note that our results for the van Hove 
singularity and the disappearance of the deformation po- 
tential are consistent with what was found in Ref. [27]. 
Although the phonon model in our work is the same, the 
dot part is different. 
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FIG. 2: (a) Spin relaxation rate as a function of magnetic field 
for the lateral length lo — 30 nm, the width a = 130 nm, and 
temperature T^lOO mK. The SO couplings Xr and Xn are 
set equal to 5 x 10~^^ and 16 x lO"'^^ eV m, respectively. The 
insets further show the enlarged regions of arrow 1 (upper 
inset) and arrow 2 (lower inset), (b) Three phonon group 
velocities vs the magnetic field, (c) The values gn and qt vs 
the magnetic field. 

Let US first focus on the dependence of the relaxation 
rates on the magnetic field B for lateral length Iq = 30 
nm. Unlike the situation in bulk system, an enhanced 
spin relaxation rate occurs as shown in Fig. 2(a) (arrow 
1 in the upper inset). This phenomenon originates from 
the van Hove singularity that corresponds to a minimum 
in the dispersion relation ujn,qi^ for finite . We further 



1.5X10''- 

|l.0x10M 
5.0x10"- 
0.0 



,200 



100 



S.O 0.5 1.0 1.5 2.0 



0.0 



0.5 1.0 
B[T] 



1.5 



2.0 



FIG. 3: Spin relaxation rate for the lateral length Iq — 60 
nm, width a = 130 nm, and temperature r=100 mK. The 
SO coupfings A_R and Ad are set equal to 5 x 10"^"^ and 16 x 
10~^^ eV m, respectively. Two enhanced and suppressed rates 
(arrow) occur. The inset shows the energy spacing AE vs the 
magnetic field B for different lateral lengths: Iq — 30 nm 
(dashed line) and Iq = 60 nm (solid line). Two horizontal 
lines in the inset indicate the corresponding energies for the 
van Hove singularity (dotted line) and the suppression of the 
rate (dashed-dotted fine). 



The relaxation rate for larger QDs exhibits a quali- 
tatively different behavior. As shown in Fig. 3, two 
van Hove singularities appear when varying the magnetic 
field. Besides, one also finds two suppressions of the re- 
laxation rate (arrow) near the singularities. We have an- 
alyzed the energy spacing between the two lowest states 
in the inset of Fig. 3. For small lateral size, the gap 
increases monotonically (dashed line). On the contrary, 
energy spacing for larger QDs shows a quite different fea- 
ture. The value initially increases as B increases. How- 
ever, after it reaches a maximum point, the energy spac- 
ing decreases with the increasing of the magnetic field B: 
although the Zeeman splitting increases with increasing 
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magnetic field, the spin-orbit interaction, on the contrary, 
tends to reduce the energy spacing between the two low- 
est levels. When the magnetic field is large enough, the 
spin-orbit effect overwhelms the Zeeman term and results 
in a decreasing tendency. Therefore, if the magnetic field 
is increased high enough, the dashed line (small QD) also 
shows similar behavior. This agrees well with the findings 
in Ref. [14]. From the inset, one recognizes that if the 
energy spacing exactly matches the specific phonon en- 
ergy (dotted line), the van Hove singularity will appear. 
For the case of a large lateral length, there are two van 
Hove singularities and two suppressions of the relaxation 
rate (dashed-dotted line). 
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FIG. 4: (Color online) Dependence of the specific energy spac- 
ings AE for the enhanced (black; mark) and suppressed (red 
mark) rates on the width a. The lateral length of the QD is 
30 nm. The Rashba constant is A_r = 5 x 10~^^ eVm and the 
Dresselhaus constant is Ad = 16 x 10~^^ eVm. 

Fig. 4 shows the specific energy spacings where rates 
are enhanced and suppressed as a function of the width. 
For the case of small widths, the enhanced rates (black 
mark) and suppressed rates (red mark) can be clearly dis- 
tinguished, and their corresponding energy spacings are 
relative large. With the increasing of the width, the en- 
ergy spacing between the enhanced and suppressed rates 
decrease monotonically. One can expect that if the width 
increases further, the system will approach the bulk sys- 
tem. This means that the van Hove singularity and the 



suppressed rate will be inhibited and eventually disap- 
pear. 

If one varies the vertical position of the dot, the rate 
will change due to different contributions from the di- 
latational and flexural waves. Accordingly, the van Hove 
singularities resulting from flexural waves will also be al- 
tered. For example, the ratio of dilatational to flexural 
wave's contribution is about 2.8:1 under the condition of 
_B = 1 T and vertical position z = 25 nm. However, if 
AE also changes, the contributions from two waves will 
also change. This is because the parameters {q\\,qi,n,<lt,n) 
of dilatational and flexural waves independently satisfy 
the dispersion relations. On the other hand, comparing 
the bulk phonons with the confined ones, the phonon- 
induced rates are roughly similar when varying the mag- 
netic field. However, there are two peculiar character- 
istics for the confined phonons. One feature is the van 
Hove singularity which results from a zero group veloc- 
ity such that an enhanced spin relaxation rate can occur. 
The second feature is a vanishing divergence of the dis- 
placement field. This will cause a suppression of spin 
relaxation rate, which is an advantage if considering the 
QD spin as a possible quantum bit candidate. 

We have studied the spin relaxation rate in a GaAs 
quantum dot embedded in a semiconductor slab, where 
an enhanced rate was found due to the phonon van Hove 
singularity. We found that at certain magnetic fields one 
enters a regime with quite the opposite characteristics, 
where a vanishing divergence of the displacement causes 
a suppression of spin relaxation rates. For larger dots 
there are multiple singularities and suppressions in the 
electron-phonon rates due to the interplay between spin- 
orbit coupling and Zeeman interaction. We believe our 
results to be useful for the understanding of spin relax- 
ation in suspended quantum dot nanostructures. Our 
findings also point at novel effects to be expected from 
future nano-scale systems where spin and mechanical de- 
grees of freedom are combined. 
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